Optical conductivity of the Prohlich polaron 
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We present accurate results for optical conductivity of the three dimensional Frohlich polaron 
in all coupling regimes. The systematic-error free diagrammatic quantum Monte Carlo method is 
employed where the Feynman graphs for the momentum-momentum correlation function in imag- 
inary time are summed up. The real- frequency optical conductivity is obtained by the analytic 
continuation with stochastic optimization. We compare numerical data with available perturbative 
and non-perturbative approaches to the optical conductivity and show that the picture of sharp 
resonances due to relaxed excited states in the strong coupling regime is "washed out" by large 
broadening of these states. As a result, the spectrum contains only a single-maximum broad peak 
with peculiar shape and a shoulder. 

PACS numbers: 71.38.-k, 02.70.Ss, 71.38.Fp, 78.30.-j 



Since the seminal work of Landau the polaron prob- 
lem, i.e. renormaUzation of quasiparticle properties due 
to coupUng to phonons, has been attracting considerable 
attention serving as a testing ground for non-perturbative 
methods. The Hamiltonian for polarons consists of terms 
describing free particle and phonon states (Planck's 
constant and electric charge are set to unity) 
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which interact through the particle-phonon coupling 

^fc-ph - E (^(q) bl 4-qak + H.c) . (2) 

k,q 

In Eqs. (|1I2|I . a^. and arc the particle and phonon 
annihilation operators in momentum space, respectively. 
The most popular models are the following two: (i) Hol- 
stein lattice polaron (characterized by the tight-binding 
dispersion law e(k) = 2t ^-^^ ^ ^ cos(fci) with hopping 
amplitude t, dispersionless optical phonons Wq = ojq, 
and short-range interaction vertex VH(q) = 7), and (ii) 
the continuous Frohlich polaron (FP) (characterized by 
e(k) = k"^ /2m, cjq = u^, and long-range interaction 

VF(q) = ii^Q \2\/2'Ka^mLOQ^ ^ ^/q; the convention is to 
measure energies in units of cjo, and length scales in units 
of l/^muo). 

Although there are many reliable methods of calcu- 
lating ground-state properties of polarons, the latter are 
hard to test directly in experiments. In contrast, the 
experimental data on optical conductivity (OC) spectra, 
(T/35(w), are readily available, but most theoretical meth- 
ods can not well treat the excited states, and, as a result, 
are less accurate in predicting (t{ljj). For decades, going 
back to early papers ^i^|4l, the polaron OC (especially 
within the Frohlich model [3,0,13,13) constantly attracts 
attention of theoretical community. 



Even in the strong-coupling limit the results on OC 
are full of contradictions. In the variational adiabatic 
treatment of polarons at strong coupling ||| one finds the 
so-called relaxed excited states (RES), i.e. quasi-stable 
states with the lattice distortion adjusted to the excited 
particle wavefunction. To contribute sharp peaks to OC 
spectra the decay rate of the lowest RES must be rela- 
tively small (otherwise the very notion of RES becomes 
vague). The decay rate for RES with the energy 0.066a^ 
fill llJ was calculated in the one-phonon approxi- 
mation Q and a sharp peak in OC was predicted for 
a > 5. Devreese, De Sitter, and Goovaerts (DSG) 01 re- 
produced this result by the expansion of the impedance 
Z{Ld) ^ l/a{Lo) within the Feynman-Hellwarth-Iddings- 
Platzman (FHIP) theory Isj. However, the validity of the 
one-phonon approximation is in doubt for the strong 
coupling limit. Also, even within the same FHIP ap- 
proximation 3], the expansion of the inverse impedance 
does not show sharp features at energies predicted by the 
impedance expansion Q. Finally, recent hig h-pr ecision 
simulations of the Lehman spectral function jl3i |. (7(w), 
did not reveal any stable excited states for FP with the 
a^-dependence of energy - peaks with weak energy depen- 
dence on a were found instead. It is worth mentioning, 
that the perturbation theory fails to reproduce data on 
g{Lo) even at very small a — 0.05. 

At present, the actual shape of OC for the Frohlich 
polaron is not known, and the crucial link connecting 
polaron theories to experiments is missing. We are not 
aware of any systematic-error free method for calculat- 
ing OC of continuous polarons. There are several non- 
perturbative approaches 0, 0, dealing with the 
properties determined by excited states. However, the 
Lanczos diagonalization technique |15| is limited to finite 
size systems, and the variational method encounters 
difhculties for long-range interactions. 

In this Letter, we show that the problem r nay be solved 
by the Diagrammatic Monte Carlo method based on 
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direct summation of Feynman diagrams for the particle 
momentum-momentum correlation function in imaginary 
time, {kfj{T)ks{0)) , with subsequent analytic continua- 
tion to real frequencies by means of stochastic optimiza- 
tion. Our technique is generic and works for arbitrary 
particle and phonon dispersion laws and interaction ver- 
tex. For the FP model we find that existing analytical 
methods severely underestimate multiphonon decay rates 
which broaden RES to the point that their contributions 
to OC overlap. Instead of sharp peaks one obtains a pic- 
ture of a broad single-maximum function with irregular 
shape and a shoulder. 

The real part of a/^si^^) is straightforwardly related to 
{kf,{r)ksm 



where (at zero temperature f3 ~ 1/T ^ oo) 



(fc^(r)fc,(0)) = 
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The solution of the integral Eq. Q is obtained by the 
stochastic optimization method of Ref. [T^l IT^ which 
allows to handle both broad and infinitely sharp fea- 
tures in the spectrum on equal footing. The (fc^(T)/c5(0)) 
correlator is diagonal in the particle momentum repre- 
sentation and is readily available (has a direct Monte 
Carlo estimator) from the statistics of Feynman dia- 
grams for the partition function, Z = Tr|e~'^^} = 

Tr {e-^««f,exp [- /^^ i/e-ph(T)dr] }, (see Fig.©. 
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FIG. 2: The low-energy part of OC for the Holstein polaron 
(circles) compared to the perturbation theory result (line); 
the high-energy tail is shown in the insert. 
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FIG. 3; Comparison between the calculated momentum- 
momentum correlator (circles) and the perturbation theory 
result (lines) for the Frohlich model. We indicate by arrows 
anomalies observed previously in the Lehman function. 



FIG. 1: A typical diagram for the polaron partition function. 
Solid (dotted) lines represent free particle (phonon) propaga- 
tors. In general, there will be diagrams with all particle prop- 
agators being accompanied by one or several phonon lines. 

Each graph generated by the expansion of the r- 
ordered exponent contributes 13~^ dT'kp(T' + T)ks(T') 
to the correlator, where k(r') is the the step-wise func- 
tion on the diagram circle. All calculations were done for 



finite but very large j3uj > 100, to ensure that exponen- 
tially small finite-temperature corrections are negligible. 
The Metropolis-type sampling of diagrams in continu- 
ous momentum-time was performed using Diagrammatic 
Monte Carlo technique developed in Refs. O 1^ ■ We 
note, that formally our method works for finite temper- 
atures as well, but analytic continuation is notoriously 
unstable at high temperatures {ijJoT ~ 1 in our case). 
To test the method efficiency and reliability for 
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FIG. 4: OC data for the weak coupling regime (open cir- 
cles) compared to the second-order perturbation theory (dot- 
ted lines) and DSG calculations (solid lines). 



the study of OC, we simulated the three dimen- 
sional Holstein polaron in the weak coupling regime 
{t = 2, ujQ = 0.1, 7 = 0.1) and compared our 
data in Fig. [21 to the perturbation theory result [l8| . 
cTaai^j) =<t(w) = {j'^/Stt'^uj) J^^ d\ qI 6{lj ~ uiq - e{q)) , 
where the integral is taken over the Brillouin zone. Not 
only the low-energy part is in perfect agreement with 
the theory, but even sharp features in the tail due to van 
Hove singularities are reproduced with remarkable (for 
analytic continuation procedures) accuracy. 

There is no such agreement with the perturbation the- 
ory for the Frohlich polaron even at extremely small cou- 
pling constant a — 0.01 (see Fig. This is due to the 
singular 1/q-dependence of the coupling constant. We 
note that positions of anomalies in a{uj) and g{uj) (see 
Ref. 0) spectra correlate. Apart from weak anomaly 
around w ~ 3, the OC spectrum is in reasonable over- 
all agreement with the second-order perturbation theory 
result for a ~ 0.01 (see Fig.^J. 

For a — 1 our results deviate from the 2-nd order per- 
turbation theory considerably; the aCTeement with the 
DSG method of impedance expansion '7] is much better. 
Since fourth-order corrections are very small in this pa- 
rameter range and do not lead to the DSG curve [lOl| . 
we conclude that two-phonon terms are not sufficient to 
account for the discrepancy and higher-order terms have 
to be taken into account even at a = 1. For a = 3 
the spectrum broadens considerably and its shape in the 
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FIG. 5: OC data for the intermediate coupling regim (open 
circles) compared to the DSG approach (solid line). Arrows 
point to the anomalies in absorption spectra arising at the 
two- and three-phonon thresholds. 



maximum starts deviating from the DSG result. 

In the intermediate coupling regime 3 < a < 7 we ob- 
serve that the position of the spectrum maximum does 
follow the DSG prediction (and that of other papers ad- 
vocating the RES concept |i[ll[Il[l3) 

for the largest 

peak, and the low-energy behavior is also captured by 
the theory reasonably well. However, the second peak at 
higher frequencies never actually fully develops, and, at 
best, its remnants are seen as a flat shoulder for a 6. 
The other qualitative difference between our results and 
DSG theory is significantly larger broadening of the dom- 
inant peak. We stress, that broadening in our data is 
not an artifact of the analytic continuation method be- 
cause we clearly see thresholds corresponding to two- and 
three-phonon emission/ absorption processes (see Fig. El . 
The conclusion we draw from this observation is that 
excited states are short-lived and their decay rates medi- 
ated b y m ulti-phonon processes are not accounted for in 
illlllllll. 

In the strong coupling regime, the broadening of ex- 
cited states completely changes the overall shape of the 
spectrum (see Fig.lSj). Instead of the double-peak struc- 
ture with extremely narrow central peak and a pseu- 
dogap, a{u!) features a broad single-maximum function 
away from the threshold. Even more surprisingly, the 
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FIG. 6; OC data for the strong coupling regime (open circles) 
compared to the DSG results (solid lines). 



shoulder at the high-frequency end of the spectrum van- 
ishes, and another shoulder emerges at the low-frequency 



end (but away from the threshold at ui = luq)- The cen- 
tral peak and the new shoulder can be fitted (somewhat 
arbitrarily) as resulting from two strongly overlapping 
Gaussian peaks. In any event, the concept of quasi-stable 
excited states clearly has no place even at large a = 8 
(the effective mass is already heavily renormalized at this 
point by a factor of thirty [131 )• 

We conclude that physically attractive picture of stable 
RES, naturally emerging from the variational treatment 
of the strong coupling limit, may not be used to interpret 
OC spectra of Frohlich polarons. Decay rates of excited 
states are severely underestimated in the one-phonon ap- 
proximation 0, . We may not rule out the possibil- 
ity that at larger values of a sharp peaks do appear in 
OC. However, this outcome would be rather academic 
because for a > 15 the effective mass renormalization is 
enormous, m*(a = 15) ~ lO'^m, and lattice effects enter 
the problem. In realistic models of polarons one also has 
to include into the picture (i) several optical modes, (ii) 
interactions with acoustic phonons, (iii) screening of the 
interaction potential, (iv) finite lattice constant, etc. The 
simulation method developed in this Letter allows one to 
study all of the above mentioned effects. Therefore the 
accurate analysis of the experimental data is possible, 
which is left for future studies. We note, that formally 
our method works for finite temperatures as well, but an- 
alytic continuation is notoriously unstable at high tem- 
peratures {ujqT ~ 1 in our case). 
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